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! I. INTRODUCTION 

a : 

. The global polarization effect in quark scatterings in non-central heavy-ion collisions has been predicted due to 

• orbital angular momenta resided in the system as a result of longitudinal flow shear P, [2, H, Q • The effect is defined 

, with respect to the direction orthogonal to the reaction plane determined by the vector of impact-parameter and 

' ■ the beam momentum. The polarization of quarks can be partially carried by hadrons containing these quarks via 

' hadronization. For example, as a consequence of this global polarization, vector mesons can have spin alignments in 

. non-central heavy-ion collisions [1] . Polarized photons are recently proposed as a good probe to the polarized quarks 

0- n . . . 

^ , In Ref. the polarization in the quark-quark scatterings has been estimated in a hot medium. The differential 
cross-section with respect to the partonic impact parameter is derived by inserting into the differential cross-section 
, in mementum space a delta function for transverse momenta, which can in turn be written as an integral over the 
partonic impact parameter. This derivation has an ambiguity. In this note we will give an alternative and a general 
way of deriving the differential cross-section with respect to the partonic impact parameter. We will show that the 
^ ■ differential cross section with respect to impact parameter in our framework is well defined and free of ambiguity. 

CN ■ II. KINEMATICS SETUP 

<N , 

The geometry of a nucleus-nucleus collision at impact parameter b is illustrated in Fig. 1 of Ref. [l|. We assume 
that two nuclei move along ±2 directions and collide at z = 0. The reaction plane is in the y direction. We consider 
the polarization along the y direction. 

Let us consider the scattering of quarks with different fiavors (Pi,Ai) -I- {P2,^2) (Pst^s) + {P4, ^4) through 
the Hard-Thermal-Loop (HTL) resummed gluon propagators, where with (i = 1,2,3,4) are spin states and Pi — 
{Ei^pizjTpiT) are 4-momenta for colliding quarks, with longitudinal and transverse momenta and p^T and energies 
Ei — pf^ + pfrp . Note that we treat all quarks massless. We assume the initial momenta are along the z direction, 
d ^ i.e. piz = —p2z > and piT — P2T — 0. We will use the shorthand notation pr = Par- The total energy in the 
center-of-mass frame is denoted by y/s = \J (Pi -I- P-i)^ . 



00 

o 



- '—I 

X 



III. AMBIGUITY OF DIFFERENTIAL CROSS SECTION WITH RESPECT TO IMPACT 
PARAMETER IN CONVENTIONAL TREATMENT 

In a conventional derivation of the differential cross section at impact parameter x-r, one starts with the one in 
momentum space, 

do = ^^\p^,p^\ E |M(PiAi,P2A2,P3A3,P4A4)P 

Ai , A2 , A4 

where a sum over the spin states Ai , A2 , A4 except A3 and an average (with a factor 1 /4) over the initial spin states 
have been taken. One can integrate out P4 by consuming the delta function 8^^^ (pi + P2 — P3 — P4) for 3-momentum 
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conservation, and then carry out the integral over p^z to remove S{Ei + E2 — E3 — E4) for the energy conservation 
which gives a factor 



J dpsJiEi +E2-E3-Ei) = ^ 



E3E4 



r,2 \E3P^:1 - Eipf, 

where p^ll — piz + P2z — and pg^j are roots of the energy conservation equation and given by 



(2) 



(1/2) 



±- 



E, 



2pl 



\PlzP2z \ -PlzP2z 



1/2 



Plz +P2z 



Then Eq. ([T]) is simphfied as 

da = 



j2 E \M{P3X3,p,\,r 



A2 , A4 



(3) 



(4) 



where we have suppressed the (PiAi,P2A2) dependence of the amphtude. Note that p4T = — Pst = — Pt is imphed 
in the above expression due to momentum conservation and the assumption that the initial state momenta are along 
the z-axis, piT = P2T — 0. One can rewrite it by inserting a delta function for transverse momenta, 



da 



1 



d^PT d'^p'rp ,^,2x(2 



64\Pi-P2\J (2^)2(27r)2 
1 



{2nf5^'\pT-p'T) 



E E I (j) (i) 
=1,2 \,.MM \E3Plz - EiP\i 



M(P3A3, P4A4)Af *(P^A3, P^A4) 



64 1 P] 



i ( ,2 d'^VT c^^Pt »(pT-p;-)-x 

\-P2\j ^{27:f (27r)2" 



E E 



M(P3A3, P4A4)Af *(P;^A3, P^A4), 



(5) 



where Pg = (-E3,P3z,Py) and P[ — (£"4,^42, — Py)- Then we obtain the differential cross section at impact parameter 

XT, 



1 



d^a _ 

d^^T ~ 64|Pi -Pil J (27r)2 (27r)2 

xE E 



i— 1,2 Ai ,A2 ,A4 
72^/^2, 



M(P3A3, P4A4)M*(P^A3, PiA4). 



(6) 



Note that the above expression of d^a /d'^'x.j' is not unique, since one can make the replacement in Eq. ([5]), for example, 

1 1 1 



E^p^ll - EJil 



E3pt - E4P); 



E^pfz - E'^P 



3z 



(7) 



with fli + 02 = 1, while keeping the total cross section unchanged. Of course one can make many other choices which 
conserve the total cross section. In Ref. ai =02 = 1/2 is used, while Eq. ^ implies ai = 1,02 = 0. So we see 
that there is an ambiguity in d'^a/d'^'x.x, or in other word, d^a /d^'x.T is not unique by this definition. The problem 
arises when one calculates quantities like the polarization when the integral over is not made in the whole space 
(therefore the delta function 5'^'^\pt — p't) is not recovered), which would lead to different or inconsistent results. 



IV. DIFFERENTIAL CROSS SECTION AT FIXED IMPACT PARAMETER 



In order to solve the ambiguity in the previous section, we will derive in this section the cross sections of parton- 
parton scatterings at fixed impact parameter in a general approach. To this end, we need to introduce particle states 
labeled by transverse positions and longitudinal momenta, which we call states in the mixed representation. They 
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are connected with states in momentum space by Fourier transform in transverse sector. We express in the mixed 
representation the final states in the scatterings as 



|P4.,A4,X4t) = / ^^^^e'P--- |P4, 



A4) , (8) 
where At is the area in the transverse plane. The S-matrix element from the initial to final states then reads 

Sfi = (P32, A3,X3t;P4z, A4,X4t| 5 |pi, Ai, P2, A2) 

^ r ATdW^4rdW^_.p3,.X3.^-.p...x.. (2^)4^(4) (p^ +P2-P,- P4) 

J (27r)^ (27r)^ 

where the final state momenta are P3 — {E^,p^z, P3t) f^nd P4 = (-£'4,^42, P4t)- The squared matrix element becomes 

2 _ ( AtcPi^^t Ard^p^T Ard^p'^T ^t^^pW -z(p3t^pU)-x3t^-»(p4t-p^^)-x4x- 
' J (27r)2 (2^)2 (27r)2 (27r)2 

; (271)8(5(4) (Pi +P2-P3- Pi)5'^^^ (Pi + P2 - P^ - Pi) 



X 



''V^ lC,E,EjE,E.EiEi ft^''"-'^3A3, PiA,), (10) 

where P3 = (-^3,^32, Pgj.) and P4 = [E'^^pj^ztV'at)- Note that the longitudinal momenta of P3 and P3 are the same, 
so are P4 and P4. The transverse parts of the delta-functions ensure P3T = — P4t = Pt and Pg^- = — P47- = p't- We 
can integrate out p^ and p-r to get rid of the four delta functions in the transverse sector and arrive at 

\Sf,\^ = A^J d2pTd'pTe^''P^"P^^'*"=^^~""^^^'°"^(Pi + P2 - P3 - P4)5'"'"^(Pi + P2 - P3 - Pi) 

^ imEjE,E.EiEi ft^')"-<^3A3. PIX.). (11) 

where ^(O-^) denote the delta functions for the energy and the z component of momenta. The differential cross section 
is then 



_ d^XT V f Ldp3z Ldpiz 

" ~ At ^VreiT ^^I2tt 2^ ' 

= ^'"^764(^ / ^^PTrfP3.d^pWP4.e-(P3-P3.) - 

X A/(-P3A3,P4A4)M*(P^A3,PiA4), (12) 

•^1 1-^2 1-^4 

where we have defined xr = X3T — X4t. In the first line we have used that the differential cross section is proportional 
to the fraction d'^xx/Ax with At = J d'^xx- Here L is the length along the z direction and Vrei = |Pi ■ P2|/(£'i£'2) 
the relative velocity of incident partons. We also used V = AtL and 2ttS'^^'^ (0) = L (because two delta functions for 
z-momenta are identical). Note that r is the time period of the scattering. It is obvious to see d^a/d^xx > from 
Eq. (HID. 

We can evaluate Eq. (|12p as follows. First we intergrate out p4z to remove d{piz +P2z —Psz —Piz)- The remaining 
two delta functions enforce energy conservation, which can be removed by carrying out integrals over the magnitudes 
of transverse momenta p'j, and px- We end up with 

rl^rr 1 1 t"^^ f'^''^ /"Pal""' P P 

' d^ d^' ^pg^g-»PT(cO.^-COS^')-T_E^4 



'Mz I r, |2 



d^xx r 64(27r)3|Pi •P2I 7o Jo Jp^r {E3 + E^)^ 

M(P3A3,P4A4)M*(P^A3,PiA4), (13) 



X 
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where and ip' are azimuthal angles of pr and relative to the direction of xt respectively. The magnitude of 
transverse momentum satisfying the energy conservation for P3 4 and P'^ 4 is proved to be the same and given by 



Vt = 



{El + E^f 



-Pi. 



[Plz +P2z -PSz)' 



n 2 



2(^1 + E2) 



1/2 



(14) 



The energies are 



E. 



E', 



pIz+Pt^ 



Ea= E'^ = y {pi,_ +P2z - PZzY + Pt- 

The integral is in the range [p™",P3^°^] where 



(15) 



PSz 



P3z 



[El + E2)^ - [piz + P2z 



(El 



2(^1- 

^E^r 



E2) 

- {Plz 



Plz + P2z 



'P2z) 



2 r 



1 + 



E1 + E2 

Plz + P2 



El + E2 



(16) 



2(^1 + E2) 
The amplitude square is given by 

Af (P3A3, P4A4)M*(P^A3, PiA4) = g\gEiE2 ^ E-iE^E'^E'^jf jf A^, (Pi - P3)A;,^,(Pi - P3O (17) 

Here g is the quark-gluon coupling constant and ctg ~ g'^/{ATi). The color factor for qq scatterings is Cqq — 

{\/N'^){5°'^5°'^ /A) = 2/9, where l/iV^ is from the average over the initial state colors. and Jj*^ are tensors 

only dependent on momentum directions, 



''2 



1 



Ei^/EiE,^ 
1 



=Tr[u(P^, A3)^I(P3, A3)7'^Pi.7"7''' 



E2\JE^E'^ A4 



= ^ Tr[u(P^, A4)m(P4, A4)7''i^2.7"7'''] 



(18) 



where u{Pi,Xi) denotes the spinor for the parton i with the spin state A; along the reference direction n = Gj,, 
and Ui,\i = u^{Pi,Xi)jo is its conjugate. A^i, is the HTL resummed gluon propagator 0, B 0- Here we only 
take the magnetic gluon exchange into account which invloves the magnetic mass /i,„ is introduced to regulate the 
divergence arising from the soft gluon exchange. The magnetic mass /i„i is proportional to the temperature T, 
Hm = 0.255{Nc/2)^/^ g^T . The numerical result from Eq. (fT5|) is given in Fig. [T]for collisions in the center-of-mass 
system of the colliding quarks. The polarized part is much less than unpolarized one. Both the unpolarized and 
polarized differential cross sections show an oscillation feature. 

Thers is an alternative way to do the integrals in Eq. (|12|) . One can intergrate out p^z to remove S{piz+p2z—P3z—Piz) 
and then p^z to remove 6{Ei + E2 — Ej, — E4) as in Eq. Wc finally obtain the differential cross section with 

respect to the impact parameter, 



J d'pTd'PTe-''''^-^'-^-'^S{E, + E,-E',-E',) 



E 



E 



1 



E3E4 



,= 1,2A„A..A. \Pl-P2WEsE,E',E', 
XA/(P3 A3, P4A4)M* (P^ A3, Pi A4), 



E3P4Z ~ E^p^^ 



Note that all final state energies are functions of px and p^. The root Pgj is given by Eq. ^ and p^^ is given 



(19) 



by pIz = Plz + P2z — Psz ■ We now rewrite the remaining delta-function into an integral over a time t which is the 
conjugate variable of the uncertainty of the final state energies arising from the specified transverse positions in the 
final state, 



2^^(P3 + P4 - S3 - k: 



•/2 



dt (,KEi+Ei-E'^-E'^)t 



T/2 



(20) 
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Figure 1: (color online) The polarized and unpolarized diflterential cross sections. The exact results from Eq. (|13|) is in the red 
solid line (unpolarized) and blue dashed line (polarized part multiplied by a factor of 10). The sum of the unpolarized part 
+ '^upoi)/d'^^T from Eqs. (|29l32l34p is in the red dotted line, while that of the polarized part dP{<yf^i + o-^^',)/d^XT 
(multiplied by a factor of 10) is in the blue dash-dotted line. The parameters are chosen to be y/s — 20 GeV, T = 200 MeV, 
r = 1.1 fm and as — 1. In the perturbation approach to obtain ^^((t'''^ + a^^^)/d^'x.T we use p^"* = 3 GeV. 
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where we have assumed T — > oo. Inserting the above into Eq. (jl9p . we obtain 



1 



E3E4 



where the time average 



)j is defined by 



1 



72 



r/2 



dt (• • •) (.KE3+Ei-E'3-E'^}t_ 



(21) 



(22) 



At small transverse momenta in the center-of-mass system where pt ^ p't ^ V^i have {E3 + E4 — E!^ — E'^) 
Pt/Vs ~ p't/V^- If Pt\^\/V^ ^ Pr'^/y/s <C 1 or r ^ V^/Pt^ we can expand the phase factor in Eq. (pij) as 



e'(E,+E^-E',^Ei)t _ ^ ^ ^(^^ + Ei^E':i- E'^)t - ^{E^ + E^ - Ef^^ - E'^ff -\ 
To the leading order, we have £'3 « E'^, E^ w E'^ and eiiE3+Ei-E'^-E'^)t ^ ^ becomes 



(23) 



1 



^ A — 1 O \ \ . 



(P^T 64(2^)4|Pi .P2I 

X Af(F3A3, P4A4)Ar(P;^A3, PiA4). 

It is interesting to see the above is unique to the leading order since 



E,pt - E,ptl EJt - E.ptl E-p't ~ E'.pf, 



Si) 



Si) 



(24) 



(25) 



The next-to-leading order non-vanishing contribution comes from the term ~ in Eq. (I23p since the linear term is 
odd in t whose integral gives zero in the range [— t/2,t/2]. It reads 



1536(2 
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xE E 

2—1,2 Ai ,A2,A4 



{E, + Ei- - E'^} 



l\2 



E'iEi^ 



\J E^EiE'^E'^ 
xM(F3A3, P4A4)Af*(P:^A3, P^A4). 



E,p^:} ~ E^p^l 



(26) 



For a simple case in the center-of-mass frame of two colliding quarks, where pi^ = —P2zi P3z = ~Piz a-nd Ei 
E2 = E3 — E4 — y/s/2. To the leading order, we have and E'^ — E'^ Ri \/s/2, then Eq. ((24|) is evaluated as 



d2xT 



367r 

X 



)-xi 



Pt + /^m P't + f^?n 



1 + iXsni ■ 



n X 



PT Pt 



9 



^^(a;T) + ^Aan • (rii x 9.t)A{xt) 



dA{xT) 
dxT 



upol , , "pol 

+ A3- 



cP 



XT 



(27) 



where we have taken only the magnetic gluon exchange into account and the magnetic mass is introduced to 
regulate the divergence arising from the soft gluon exchange. We denoted rii — as the direction of pi and used 



r 1 rf 

A{xt) = ^ = 2-/ 

J Pt ^ Mm Jo 



PtJo[ptXt) 
dpT — — 5 — : 
Pt + Aim 



(28) 



where Jq is the Bessel function and p^* is the cutoff to regulate the ultraviolet divergence. The polarized and 
unpolarized differential cross sections can be obtained 



Note that the above result is just Eqs. (40,41) in Ref. 
The next-to-leading order contribution is evaluated as 



2"'„.(„,xxt)A(xt)-'^^("^) 



dxi 



(29) 



2167r2s 
1 



1 + iX^ni 



1 



Pt + Mm p'% + 
(1) 



n X 



Pt - Pt 



d2xT 



+ A; 



d2xT ' 



where we have used 



(£;3+^4-^3-^4)' = 



n 2 



V^-2Js/4-p2,+p'2^) 



(p| - p't 
s/4 



2^2 



(30) 



(31) 



The unpolarized part reads 

d2a(i) 



upol 



d'^TiT 



216 



2 „'2 s 
T) 



\ f d2pTd2p^e-*(P--P-)-^-(pf, + p'^T - ^PtP' 
n^s J 

(AoAi - Aj), 



1 



1 



Pt + Mm Pt + M^ 



r2a| 
27s 



where we used 



A,{xt) 



Pt 



dpi 



Pt 
PT + y-l 



■Jo{ptXt) 



(32) 



(33) 
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Figure 2: The polarized and unpolarized differential cross sections in the leading order from Eq. (|29|) (left panel) and the 
next-to-leading order from Eqs. p2l34|l (right panel). The parameters are set to the same values as in Fig. [1] 



0.20 

0.15 
0.10 I- 
0.05 
0.00 



0.2 



Unpolarized 
Polarized 



0.4 0.6 
xt (fm) 



0.8 1.0 



0.01 
0.00 
^^-0.01 
-0.02 
-0.03 





/ 

/ 
/ 

1 

1 

1 

1 


Unpolarized - 


1 


Polarized 







0.2 



0.4 0.6 
Xt (fm) 



0.8 1.0 



For i = 0, 1, 2 with rii — 1, 5, 3. The polarized part turns out to be 

The numerical results of Eqs. (|29l32l34p are shown in Fig. [51 We see that both the leading and next-to-leading parts 
are damped out above 0.4 fm. The next-to-leading contributions of the unpolarized and polarized parts are about 
1/6 of the leading counterparts. Both the unpolarized and polarized differential cross sections show an oscillation 
feature. The sums of the leading and next-to-leading contributions are shown as the red dotted (unpolarized) and blue 
dash-dotted line (polarized) in Fig. [TJ We already mentioned that a cutoff in transverse momentum "* is needed to 
regulate the integrals in the leading and next-to-leading differential cross sections. The cross section results depend 
on in the perturbation. The time scale r is also a quantity to be determined. We can find the range of the r by 
equating the exact result d^a/d^XT from Eq. and the sum <P{a^°^ + a^^'')/cPxT in the perturbation approach 
from Eqs. (|29l32l34p at a specified value of p^"*. 



V. SUMMARY AND DISCUSSION 



We have proposed a general approach to the differential cross section with respect to impact parameters, which is 
well-defined and free of ambiguity existing in the conventional approach. The main difference of our approach from 
the conventional one is that (1) in the conventional approach the transfer of small transverse momenta is implied, 
while the general approach is valid for all transverse momenta; (2) there are two independent delta functions in the 
general approach for energy conservation in the cross section formula, which arises from fixing impact parameters in 
the final state partons making the total final state energy be uncertain. While in the conventional approach the two 
delta functions for energy conservation are identical turning the second one to be the infinite interaction time t. 

As a simple illustration of our formalism, we evaluated the polarized and unpolarized differential cross sections at 
small angle quark-quark scatterings in the center-of- mass system of the colliding quarks. To smoothly connect the 
general approach and the conventional one, we propose an expansion in terms of AE — Ef — Ei ~ l/r with Ei and 
Ef the initial and final state energies in collisions. The leading order contribution reproduced the conventional result, 
i.e. that of Ref. Q. 

The general formulation in this paper can also be applied to many other parton-parton scatterings in heavy ion 
collisions or even proton-proton collisions [Toj . 
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